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Abstract 


Background: Migella sativa, generally known as black seed or black cumin, has a rich 
historical background in traditional medicine for its diverse health benefits. Objective: The aim 
of the study explore the potential effects of Nigella sativa (NS) on colorectal cancer (CC) 
through the application of network pharmacology and molecular docking. Methods: Network 
pharmacology (NP) analysis was conducted to identify pertinent colorectal cancer targets and 
compounds sourced from relevant databases. Subsequently, protein-protein interaction (PPI) 
and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were employed to 
delineate critical molecular pathways, while molecular docking simulations study were 
employed to validate the binding interactions at active sites. Results: Network analysis of 
compound-target interactions revealed a network comprising 319 nodes and 410 edges, 
indicating a complex interplay between the derivatives and their associated targets. The PPI 
analysis underscored the significant interactions within the network, particularly with targets 
known to be involved in colorectal cancer regulation. KEGG pathway analysis highlighted the 
importance of EGFR and PI3K pathways in the context of colorectal cancer. Notably, molecular 
docking studies identified Nigellicine as having the highest affinity for key colorectal cancer- 
related targets, including AKT1 (-7.6 kcal/mol), IL6 (-6.2 kcal/mol), ALB (-7.1 kcal/mol) and 
HSP90AA1 (-8.2). Conclusion: The integration of network analysis and molecular docking 
studies provided collectively support the colorectal cancer characteristics of the compounds, 
emphasizing the need for further research to develop novel pharmacological interventions for 
colorectal cancer management. 
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Introduction colorectal cancer ranks third in terms of prevalence 


worldwide [2]. Hereditary disorders such as familial 


Cancer is characterized by excessive cell proliferation that 
can metastasis to other parts of the body and cause disease. 
The rising population, the aging population, and the effects 
of lifestyle changes all contribute to a rise in cancer 
incidence and mortality rates. Because cancer has such far- 
reaching effects, it degrades the quality of human resources 
[1]. As the second biggest cause of cancer-related death, 


adenomatous polyposis (FAP) and hereditary nonpolyposis 
colorectal cancer (HNPCC) account for 10% of all cases of 
colorectal cancer [3]. The remaining 70% of cases arise at 
random. 

Colorectal cancer treatment is complex and calls for input 
from other fields of expertise. There is a higher recurrence 
rate with local operations, surgical therapy is less successful, 
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especially for metastatic disease, and the adverse effects of 
chemotherapy and monoclonal antibodies are distinct [4]. 
This issue is mostly attributable to the fact that the 
medications being utilized are not selective, causing harm to 
both cancer cells and healthy cells. Numerous sets of genes 
are involved in most disorders. However, creating a new 
medication is a lengthy and expensive process [5]. Because 
of this issue, the focus of drug development has shifted from 
developing single-drug therapies to creating plant-based 
medicines that utilize many chemical compounds directed at 
different targets (multicomponent - network targets) [6]. 
Multicomponent medications have more than one active 
ingredient, and this can have a positive impact on the 
pharmacological effects. In silico computational methods 
have a synergistic effect on multicomponent approaches, 
allowing for faster execution and more encouraging 
outcomes. The concepts of network theory between 
chemicals and biological systems, as well as the outcomes of 
in vitro and in vivo testing for drug development, are 
explained by network pharmacology [7]. 

For correct findings and productive preclinical research, an 
early investigation providing an initial view is crucial. 
Compound content and mechanism of action are predicted 
using computational methods based on the compound's link 
to proteins in existing databases. The purpose of this 
research is to foretell Nigella sativa network pharmacology 
in colorectal cancer. Our findings can be used as a 
foundation for developing safe and effective new medicines 
or improvements in cancer treatment based on target genes. 


Materials and Methods 


Screening for active components of Nigella sativa 

The constituents employed in this investigation were 
procured from the KNApSAcK Family database 
(http://www.KNApSAcKfamily.com/) and Dr. Duke’s 
Phytochemical and Ethnobotanical Databases 
(https://phytochem.nal.usda.gov/phytochem/search). The 
selection of active compounds adhered to specific criteria, 
including a molecular weight (MW) within the range of 180 
to 500 Dalton, oral bioavailability (OB) equal to or greater 
than 20%, drug-likeness (DL) equal to or greater than 0.1, 
and a blood-brain barrier (BBB) value equal to or greater 
than —0.3. The chemical structure validation of Nigella 
sativa components was conducted utilizing the PubChem 
database (https://pubchem.ncbi.nlm.nih.gov) and was 
generated using ChemDraw. 


Screening for protein targets of Nigella sativa 

The protein targets associated with Nigella sativa 
constituents were sourced from the Swiss Target Prediction 
database — (http;//www.swisstargetprediction.ch/). Target 
information relevant to colorectal cancer was extracted from 
the Gene Card Database (https://www.genecards.org/) and 
the Online Mendelian Inheritance in Man (OMIM) database 
(https://www.omim.org/). Putative targets of the Nigella 


sativa compounds were identified through a screening 
process that involved detecting common targets among the 
genes represented in a Venn plot 
(https://bioinformatics.psb.ugent.be/webtools/Venn/). 


Protein-protein interaction (PPI) data 

The exploration of protein-protein interactions (PPI) serves 
as a valuable approach to discern potential therapeutic 
targets in the management of disorders, given the pivotal 
role of PPI in orchestrating biological processes [8]. In this 
investigation, the shared targets of Nigella sativa compounds 
and colorectal cancer were delineated utilizing the Search 
Tool for the Retrieval of Interacting Genes/Proteins 
(STRING) V11.5 platform (https://string-db.org/) designed 
for multiple proteins in Homo sapiens [9]. A stringent 
minimal interaction score of 0.9 (indicative of high 
confidence) was applied in the establishment of the PPI 
network. Subsequent extraction of p-values was conducted 
to assess the statistical significance of the PPI network and 
its enrichment. 


Gene and KEGG enrichment analysis 

The investigation of signal pathways, biological processes, 
molecular functions, and cellular component terms 
associated with the intersecting targets between the drug and 
disease was conducted through Gene Ontology (GO) 
functional annotations and the Kyoto Encyclopedia of Genes 
and Genomes (KEGG), considering a significance threshold 
of p « 0.05. Enrichment analysis of these pathways was 
performed using Input 2.0 
(http://cbcb.cdutcm.edu.cn/INPUT/Home/) and  Shinygo 
(http://bioinformatics.sdstate.edu/go/) as analytical tools [10- 
11]. 


Compound-target network construction 

The compounds derived from Mgella sativa were 
systematically linked to construct а compound-target 
network. Cytoscape 3.9.1 (Cytoscape Consortium, San 
Diego, CA, USA) [12], a software tool designed for 
visualizing interaction networks, was employed for the 
network construction. In this network, edges denote the 
interactions between the compounds and their corresponding 
protein targets, while nodes represent the individual 
compounds and their associated targets. 


Molecular docking analysis 

AutoDock Vina was employed to assess the interactions 
between Nigella sativa compounds and target proteins [13]. 
The three-dimensional (3D) structures of the target proteins, 
in Protein Data Bank (PDB) format, were retrieved from 
(https://www.rcsb.org) for subsequent docking analysis [14]. 
Additionally, the three-dimensional (3D) structures of 
Nigella sativa compounds were sourced from the PubChem 
database (https://pubchem.ncbi.nlm.nih.gov). AutoDock and 
PyMOL software were used to automatically determine 
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docking scores (binding affinities) and binding interaction 
between the compounds and target proteins [13]. 


Results 


Screening for active components of Nigella sativa 

Ten components inherent in Nigella sativa were anticipated 
through a comprehensive analysis utilizing the KNApSAcK 
Family and Dr. Duke's Phytochemical and Ethnobotanical 
databases. Specifically, five components were identified 
within the seed oil of Nigella sativa L., as detailed in Table 
1. The KNApSAcK database, referenced in [15], provides a 
wealth of information encompassing the interrelation of 
components, biological activities, and associated species. 
Additionally, Dr. Duke's Phytochemical and Ethnobotanical 
database, accessible online at no cost, offers valuable ethno 
medicinal data of high significance. Nigella sativa seeds 
have historically played a pivotal role in treating various 
diseases, including cancer, within traditional medical 
systems prevalent in South and Southeast Asia, Arabia, 
Africa, and the Mediterranean [16]. Employing 
pharmacokinetic criteria and literature reviews, we identified 
nine potentially active chemicals from Nigella sativa, 
detailed in Table 2. These compounds exhibit molecular 
weights (MW) ranging from 180 to 500, oral bioavailability 
(OB) between 20 and 30, and blood-brain barrier (BBB) 
binding affinity constants between 0.3 and 0. The chemical 
structures of these selected compounds were delineated 
using ChemDraw, as illustrated in Figure 1, showcasing a 
diversity of structures encompassing glucosides, long-chain 
lipids, and simple chemical molecules. 

The study reveals the capacity of Nigella sativa components 
to interact with a diverse array of proteins, including those 
implicated in cancer, as determined through protein target 
prediction utilizing the Swiss Target Prediction database. 
The database identified 410 potential chemical targets 
associated with Nigella sativa. To identify potential targets 
pertinent to colorectal cancer, the Gene Card Database was 
employed. The combined analysis of Swiss Target 
Prediction and GeneCards Database yielded а 
comprehensive list of Nigella sativa compound protein 
targets, as detailed in Table 3. 

The selection of 410 targets predicted by the Swiss system 
was visually represented using Cytoscape 3.9.1, resulting in 
a compound-target network displayed in Figure 2A-2B. The 
Venn diagram illustrates shared targets and genes, involving 
9 chemicals and 410 interacting target proteins. Notably, 
multiple components within this network were observed to 
target the same proteins. The network comprises 319 nodes 
and 410 edges, suggesting a potential synergistic effect of 
Nigella sativas active biochemicals on various targets. This 
aligns with the plant's documented therapeutic efficacy 
against a spectrum of diseases and ailments. Furthermore, 
topological characteristics, including betweenness centrality, 
closeness centrality, and degree, were assessed to elucidate 
the importance of nodes within the network, as summarized 


in Table 4. These findings underscore the multifaceted 
potential of Nigella sativa compounds in targeting specific 
proteins, thereby presenting a promising avenue for 
therapeutic interventions against diverse diseases, including 
colorectal cancer. 


PPI network 

Figure 3A-B illustrates the protein-protein interaction (PPI) 
network depicting potential colorectal cancer targets 
associated with Nigella sativa. The PPI enrichment analysis 
yielded a p-value of 0.0272, while the average node degree 
was determined to be 8.6. Notably, specific targets such as 
АКТІ exhibited a node degree of 120, IL6 demonstrated a 
node degree of 112, ALB exhibited a node degree of 110, 
and HSP90AA1 showed a node degree of 88. These 
findings suggest that these identified targets, particularly 
АКТІ, IL6, ALB, and HSP90AA1, may play pivotal roles 
in mediating the effects of Nigella sativa on colorectal 
cancer. 


GO and KEGG enrichment analysis 

Genes with potential relevance to colorectal cancer 
treatment were identified through a Gene Ontology (GO) 
enrichment study of Nigella sativa. These genes were 
associated with various cellular components, molecular 
functions, and biological processes. Noteworthy biological 
processes included positive regulation of gene expression, 
response to endogenous stimulus, cellular localization, 
intracellular signal transduction, nitric oxide biosynthetic 
process, protein phosphorylation, epidermal growth factor 
response, apoptotic process, and acute inflammatory 
response, as indicated in Table 5. Molecular functions 
associated with these genes encompassed nitric-oxide 
synthase regulator activity, enzyme binding, binding of 
identical proteins, binding of phosphatases, ions, proteins, 
protein phosphatases, and small molecules, detailed in Table 
6a. Additionally, Table 6b highlighted the predominant 
cellular component term as the Endoplasmic reticulum 
lumen in the GO analysis. 

Furthermore, the Kyoto Encyclopedia of Genes and 
Genomes (KEGG) analysis Table 6c and Figure 4 KEGG 
Pathways revealed several pathways linked to the potential 
colorectal cancer impact of Nigella sativa. These pathways 
included Pathways in cancer EGFR tyrosine kinase 
inhibitor resistance, the IL-17 signaling pathway, prostate 
cancer, the ErbB signaling pathway, the TNF signaling 
pathway, the Estrogen signaling pathway, gastric cancer, 
Th17 cell differentiation, the MAPK signaling pathway, 
apoptosis, and the Rap! signaling pathway (Figure 5A-E). 
The implicated genes within these pathways encompassed 
the epidermal growth factor receptor, nitrogen metabolism, 
cancerous pathway, PI3K, and vascular endothelial growth 
factor. 
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Figure 1. Structures of Phytoconstituents of Nigella sativa L. 


Figure 2A. Compound- Target Network. 
Targets 


Genes 


Figure 2B. Venn diagram shows the intersection between Nigella Sativa targets and colorectal cancer related 
genes. 
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Figure 3. A) PPI network of potential colorectal cancer targets of Nigella Sativa. The PPI enrichment p-value 
was 0.0272 and the average node degree was 8.6. The node degrees of AKT1 was 120, IL6 was 112, ALB was 
110, and HSP90AA1 88 was suggesting that these targets could play important roles in the effects of Nigella 
Sativa on colorectal cancer. B) Top 10 Hub Genes Network. 
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Figure 4. KEGG Plot. 
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Figure 5A. EGFR Pathway. 
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Figure 5B. Nitrogen Metabolism. 
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Figure 5D. VEGF. 


16 


Ansari Vikhar Danish Ahmad et al., JIPBS, Vol 11 (1): 10-24, 2024 


Ful TED Carinogers 
PATHWAYS IN CANCER. Caen ute? Dots a ron s vA ] 9 
PITCH! kl so Notch Fas s TOPR 
Oxidative or a T { Ж 
Hedgehog signaling Notch signaling ШТ: momen + | 
vc Sie Fumante НЕН | Emgac|CUI2] "woe zt. Nop (миъ 1 Gano | fm 
* FH VHL | Rbxt Q 9.0, CASP9 D | Enos ын 
Microtubule +— KiF; [5080 AHNE On CASPR 
= Qo At) Y | TGF-B signaling 
Pam " » == Сус ]<— Mitochondrion ««—[ pid | таву GSTPL 
= 3 | ` 
E „ „ Bak 
DCC РР] Se | S МОХА Ј[РОМА ] | I 
` PEN ! Bela lx 
м | Ва m | 
Т | 
| EVIL 
ja-cateran | BRCA2 Man | [АМ SOU 
ECAD | Rad51 ҺМ5Н2| | 
1 
| DNA DNA | СВР 
! О. | Y 
e ES | | 4 | 
x * | Вах ][һМ5Н3] boeenl ^ | | 
Мєл nk | LE Lu H 
T6) GER. |в ions OSTs | [ret / E 
| 1 / — DNA damage Кос Schu 
ze >l orcr tv «i L 1 4 "d | pi 
^ 3 signal 
Q—*| / "e | RN 
PL o >| ei d Genomic instability $ 
PGE2 popes МПЕ]———> ——p| Resistance to 
Z if 3 chemotherapy 
N f (d | vá 
TIGA / / pos] ^ ач 
км ]— "roe \ /. Jj “ 
NN / / Be ][ Bak А 
"Mareceplor ( Focal sesion ) к —“ Гыз [Fak 
— er ey ы | 
ad К 
S B i Ко | 7 РА SC д Telomeric 
= RAR] Ç / P pg DNA repeats 
2 CRK " P | j m 
BK рэкан NW ва [а] 2 >x А] Poy ар [сив | 
= En Signaling pathway uet P гу, Á À 7 / E [eG Hm Bi 
rm Zei 3 Z Ji! Ween 1 .- / Z VIA if Ss m 
ai KM D OO 
IP AG E + 4 d N » TRN 
SI ~- | _____ 
N кай palway / = Viz) 
Cytokine-cytokine "WE Po 
VH, > [STATS mE Lo } Ë SCF 
mem Deag ud I Se 
> "1 stati ]J "HE / k i 
TGFs EGFR fF NC 
ЕЕЕ 2329 2 ` 
ec Ms]-—----------- p 
ЕСЕ [17:3] bag E Agi] 244 d E | [GcseR][McsFr] EN [ve] [тк 
Sag ie MECH М ВАВА 
roor —| Kops — zr |! Ex) S. tnget genes Шш genes быш: Барые 
is | бк N „к= E | SS `. t t 1 
mmc} | — cx: tel Se MEK >i < [us = [omm 1 
[zrscc] — El ll E Tz tt РА А 
pena K ; 01 
na — EH - | MAPK аныш | г Ei Esc) I 
Fá унчу = = CDKA iL put] [rui] [сева] [сева AML 
SS ]—[ FGFR. ^^ Eed мка SS = 
Fa ` с-Мус N 
f 
Z Cdim |/ \ 
„7 signaling pathway | у E Е Вавр |— Сова — HUE 
TRK kŠ 
e O| Ek] À щш) d l 
mo, K DAG \ РА (ege Tr jH 
F 
= `- somiem] Mahn Fai Ó o 
[ramis 5 ux Ligands — Fo acid 
Data on KEGG graph 
Rendered by Pathview 


Figure 5E. Pathways in Cancer. 


Table 1. Prediction of components contents of Nigella Sativa with KNApSAcK Family and Dr. Duke’s 
Phytochemical and Ethnobotanical Databases. 


Compound Name CAS ID SMILES Plant Part 

m-Thymol 89-83-8 с1сс(с(сс1С)О)С(С)С Seed Oil 

Oleic acid 112-80-1 CCCCCCCC/C-CXCCCCCCCC(-O)O Seed 

Myristicin 607-91-0 c12c(e(cc(c1)CC-C)OC)OCO2 Seed Oil 

p-Cymene 99-87-6 с1сс(ссс1С)С(С)С Seed 

beta-Amyrin 559-70-6 CC1(CCC2(CCC3(C(=CCC4C3(CCC5C4(CCC(C5( Seed Oil 
C)C)O)C)C)C2C1)C)C)C 

Linolenic acid 463-40-1 OC(CCCCCCC/C=CC/C=CC/C=CCC)=O Seed 

Thymoquinone 490-91-5 C1=C(C(=O)C=C(C1=O)C)C(C)C Seed Oil 

Kaempferol 197250-98-9 C1=CC(=CC=C1C2=C(C(=O)C3=C(C=C(C=C302) Seed Oil 

3-glucosyl-(1->2)- O)O)OC4C(C(C(C(O4)CO)O)O)OCSC(C(C(C(O5)C 

galactosyl-(1->2)-glu O)O)O)OC6C(C(C(C(O6)CO)O)O)O)O 

coside 

Quercetin 197250-97-8 C1=CC(=C(C=C1C2=C(C(=O)C3=C(C=C(C=C3O2 Seed 

3-glucosyl-(1->2)- )O)O)OC4C(C(C(C(O4)CO)O)O)OC5C(C(C(C(O5) 

galactosyl-(1->2)-glu CO)O)O)OC6C(C(C(C(O6)CO)O)O)O)O)O 

coside 

Quercetin 197294-29-4 COC1=C(C=CC(=C1)C=CC(=O)OCC2C(C(C(C(O2 Seed 

3-(6””-feruloylgluco )OC3C(C(C(OC3OC4C(C(C(OC4OC5=C(OC6=CC( 

syl)-(1->2) =CC(=C6C5=0)0)O)C7=CC(=C(C=C7)O)O)CO)O) 
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-galactosyl-(1->2)-gl O)CO)O)O)O)O)O)O 
ucoside 
Nigellicine 98063-20-8 CC1=CC(=O)C2=C(N3CCCCN3C2=C1)C(=O)O Seed 
Nigellidine 120993-86-4 CC1=CC(=O)C2=C(N3CCCCN3C2=C1)C4=CC=C( Seed 
C=C4)O 
Nigellimine 4594-02-9 CC1=NC=CC2=CC(=C(C=C12)OC)OC Seed 
Carvone 99-49-0 CC1=CCC(CC1=O0)C(=C)C Seed 
Nigellidine 1032622-86-8 CC1=CC2=C(C(=C1)OS(=O)(=O)O)C(=[N+]3N2C ` Seed 
4-O-sulfite CCC3)C4=CC=C(C=C4)O 
Table 2. Active components of Nigella Sativa. 
Molecule Name MW (g/mol) B (%) BBB DL 
m-Thymol 150.22 0.55 Yes 0.95 
Oleic acid 282.46 0.85 No 1.0 
Myristicin 192.21 0.55 Yes 0.65 
p-Cymene 134.22 0.55 Yes 0.94 
beta-Amyrin 426.72 0.55 No 0.98 
Linolenic acid 278.43 0.85 Yes 1.0 
Thymoquinone 164.20 0.55 Yes 0.93 
Kaempferol 3-glucosyl-(1->2)- 772.66 0.17 No 1.0 
galactosyl-(1->2)-glucoside 
Quercetin 3-glucosyl-(1->2)- 788.66 0.17 No 1.0 
galactosyl-(1->2)-glucoside 
Quercetin 964.83 0.17 No 1.0 


3-(6””-feruloylglucosyl)-(1->2) 
-galactosyl-(1->2)-glucoside 


Nigellicine 246.26 0.85 Yes 0.77 
Nigellidine 294.35 0.55 Yes 0.88 
Nigellimine 203.24 0.55 Yes 0.89 
Carvone 150.22 0.55 Yes 0.63 
Nigellidine 4-O-sulfite 375.42 0.55 No 0.53 

Table 3. Target proteins of Nigella sativa L. with Swiss Target Prediction. 
Compound name PubChem ID Target genes linked to cancer 
m-Thymol 6989 TRPA1, PTGSI, GABRAI, GABRG2, HTR2B, GABRB3, 


HTR2C, CA2, CHRM2, FLT3, ЈАКІ, JAK2, PRKCA, AURKA, 
ESRRG, CDK2, CCNAI, CCNA2, ACHE 

Myristicin 4276 CDKSR1, CDK5, CDK2, ССМАІ, CCNA2, CDK9, CCNTI, 
DYRKIA, МЕТАР2, PDE7A, АРОКА2В, ADORA3, 
ADORA2A, DYRK1B, GRM4, JAKI, JAK2, CHRNA3, PDESA, 
JAK3, NUDTI, AKTI, MAPK3K14, GRM5, TNKS2, NAMPT, 
ERBB2, СМКІ, HSP90AA1, ALPL, HTR2A, HTR2C, LIMKI, 
CDK2, CDC7, GABRB3, GABRA3, GABRG2, GABRAI, 
GABRAS, GABRA2, CCNEI, CDK2, CDK5, CHEKI, KAT2B, 
KDR, HMOXI, PTPRC, STS, NOSI, NOS2, NOS3, TERT, 
CD38, TGM2, TYMP, CDK1, TSPO, MAPK8, GSK3B, MET, 


MAPK10 
p-Cymene 7463 CYP2A6, ACHE, TAARI, PPARA, PTGS1, TRPAI 
Linolenic acid 5280934 PPARG, PPARA, PPARD, FABP4, FFARI, FABP3, PTGSI, 


FABPS, TERT, ЕАВРІ, CNR1, PTPNI, ALOXS, PTPN2, 
PTGES, LTB4R, POLB, ESR2, PTPN6, RORC, TOP1, PTGER2, 
ALOXI2, CDC25A, PTPRF, PTGS2, NOS2, PGR, SRD5A2, 
PDE4D, PTGERI, PLA2GIB, ЕЅКІ, CDC25B, NR3CI, 
PTPN11, ADORA3, NRIH3, CD81, PRKCH, FNTA, FNTB, 
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Thymoquinone 10281 
Nigellicine 11402337 
Nigellimine 20725 
Carvone 7439 


RXRB, SERPINA6, SHBG, G6PD, ITGAL, ICAMI, ITGB2, 
PTGER4, CYP26A1, NPCIL1, CYP17A1, AKR1B10, IMPDH2, 
1L6, HTR2B, MAPK14, MAPKI, PREP, AR, MDM2, FFAR4, 
ENPP2, MAPK3, RBP4, PRKAGI, PRKABI, PRKAA2, 
ALOXI5, TOP2A, КОКА, MMP2 

PLK1, GLD, GLI1, ALOXS, CYP19A1, PTPN2, CHRM2, CA2, 
ACHE, SHBG 

HTR2A, ADRAID, HTR2B, ADRAIA, HTR2C, ADRAIB, 
TSPO, JAK3, DBF4, CDC7, MMP3, HPGD, DYRKIA, 
CDKS5RI, CDK5, AR, EGFR, KDR, AURKA, АКТІ, SIGMARI, 
CSFIR, MMP13, MMP2, MIF, ALOXI5, ALOX12, СМКІ, 
MMP9, FGFRI, DYRKIB, GABRA1, CCNE2, CDK2, CCNEI, 
CDK1, CCNBI, CCNB2, PIMI, MAPKAPK2, GSK3B, CTSK, 
CDK2, CDK1, GABRAS5, ROCKI, ЕКМІ, AKRIBI, JAK2, 
CA12, AKRIBIO, RBP4, PIK3CG, PIK3CA,TGFBRI, SIRT2, 
RORC, MAPKI, CASP3 

ACHE, NQOI, CYP1A2, GRM4, CCNE2, CDK2, CCNEI, 
DYRKIA, DYRKIB, JAKI, JAK2, CTSK, GSK3B, CSNKID, 
GRM5, ADORA2B, IDO1, PDE5A, CDK1, CCNB1, CCNB2 
CYP19A1, SRD5AI, NR3CI, PGR, SERPINA6, SHBG, FABPI, 
NRII3, SRD5A2, AR, ADHIA, SIGMARI, ADHIC, NPCILI, 
PTGES, CA2, CYP17A1, PARPI, ADORA3, MAPK3, PRKCH, 
PTPN11, AKRIBIO, PTPN2, PTGSI, NR1H3, TOP1 


Table 4. Topological analysis of important nodes with network analyzer results. 


Name Closeness Betweenness Degree Average Shortest Topological Coefficient 
Centrality Centrality Path Length 

MAPK3 0.818182 0 7 1.2222222222222 1 

ALB 1 0.007937 9 1 0.950617 
IL6 1 0.007937 9 1 0.950617 
АКТ1 1 0.007937 9 1 0.950617 
GSK3B 0.9 0 8 1.11111111111111 0.986111 
PTGS2 1 0.007937 9 1 0.986111 
ESRI 1 0.007937 9 1 0.950617 
PPARG 1 0.007937 9 1 0.950617 
HSP90AA1 1 0.007937 9 1 0.950617 
ЕКВВ2 0.9 0 8 1.11111111111111 0.986111 


Table 5. Go Biological Process. 


Description Count in Gene Set False Discovery Rate 
Positive regulation of gene expression 9 1.18E-06 

Response to endogenous stimulus 9 2.75E-06 

Regulation of cellular localization 8 7.27Е-06 

Regulation of Intracellular signal transduction 9 1.13E-05 

Positive regulation of nitric oxide biosynthetic process 4 1.85E-05 

Positive regulation of protein phosphorylation gi 2.84E-05 

Response to epidermal growth factor 3 0.00031 

Regulation of acute inflammatory response 3 0.00035 
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Table 6a. Go Molecular Function. 


Description Count in Gene Set False Discovery Rate 
Nitric-oxide synthase regulator activity 3 7.64E-05 

Enzyme binding 8 0.0015 

Identical protein binding 8 0.0015 

Phosphatase binding 4 0.0026 

Anion binding 7 0.0343 

Protein binding 10 0.0372 

Protein phosphatase binding 3 0.0372 

Small molecule binding 7 0.0372 


Table 6b. Go Cellular Components. 


Description 


Count in Gene Set False Discovery Rate 


Endoplasmic reticulum lumen 


4 0.0258 


Table 6c. KEGG Pathways. 


Description C 


ount in Gene Set False Discovery Rate 


Pathways in cancer 9 
ЕСЕК tyrosine kinase inhibitor resistance 
IL-17 signaling pathway 

Prostate cancer 

ErbB signaling pathway 

TNF signaling pathway 

Estrogen signaling pathway 

Gastric cancer 

Th17 cell differentiation 

MAPK signaling pathway 

Apoptosis 


5 
5 
5 
4 
4 
4 
4 
3 
3 
2 
Rapl signaling pathway 2 


2.04E-11 
4.59E-08 
6.83E-08 
6.97E-08 
2.03E-06 
5.04E-06 
8.48E-06 
1.10E-05 
0.00015 
0.0019 
0.0074 
0.0147 


Molecular docking analysis 

The study conducted molecular docking analyses to assess 
the binding affinities of active constituents from Nigella 
sativa with associated protein targets in the protein-protein 
interaction (PPI) network, including AKT1, IL6, ALB, and 
HSP90AAI. Molecular docking scores, indicative of the 
strength and stability of chemical-protein binding, were 
employed, where a more negative score indicated a stronger 
interaction. Table 7 presents the docking scores for the 
interactions between Nigella sativa compounds and their 
respective protein targets. Figures 6A-D illustrate diverse 
binding mechanisms between the compounds and proteins. 
Crystal structures from the Protein Data Bank (PDB) for 
АКТІ (PDB: 4EKL), 116 (PDB: 4757), ALB (PDB: 
6МАВ), and HSP90AAT1 (PDB: 3001 were obtained. Auto 


Dock Vina was employed for docking ligands into receptors 
to anticipate potential binding interactions, with the grid 
centered on amino acid residues surrounding the active sites. 
Free binding energies, expressed in kcal/mol, were utilized 
to convey the results. Molecular docking experiments 
revealed estimated free energy of binding ranging from -5.7 
to -7.6 kcal/mol for АКТІ (PDB: 4EKL), -4.3 to -6.2 
kcal/mol for IL6 (PDB: 4ZS7), -5.4 to -7.1 kcal/mol for 
ALB (PDB: 6M4R), and -7 to -8.2 kcal/mol for HSP90AA1 
(PDB: 3001). Nigellicine emerged as the most active 
molecule, exhibiting a binding energy of -8.2 kcal/mol for 
HSP90AA1 (PDB: 3001), -7.6 kcal/mol for АКТІ (PDB: 
4EKL), -6.2 kcal/mol for IL6 (PDB: 4757), and -7.1 
kcal/mol for ALB (PDB: 6MA4R). 
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4ekl 


3001 


PHE-74 


CYS-73 


GLN-183 


4257 
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6M4R 


Figure 6A-D. Binding modes of ligands, and receptors by molecular docking simulation. Docking pose and 
surface representation of Nigellicine, and AKT1 (PDB: 4EKL). Docking pose and surface representation of 
Nigellicine, and HSP90AA1 (PDB: 3001). Docking pose and surface representation of standard drug 
Nigellicine, and IL6 (PDB: 4ZS7) and Docking Pose and surface representation of ALB (PDB: 6M4R). 


Table 7. Docking scores of Nigella Sativa compounds and potential targets. 


Binding Energy (kcal/mol) 


Compound АКТ1 IL6 ALB HSP90AA1 
(PDB:4EKL) (PDB:4ZS7) (PDB:6M4R) (PDB: 3O0D 

Nigellicine -7.6 -6.2 -7.1 -8.2 
Nigellimine -6.5 -5.3 -6.7 -7.9 
Thymoquinone -6.5 -5.1 -6.7 -7.5 
Myristicin -6 -5.2 -6.2 -7.4 

Cymene -5.9 -4.9 -6 -7.3 

Thymol -6.2 -5 -6 -7.3 

Carvone -6.1 -5 -5.8 -7.2 

Linolenic acid -5.7 -4.3 -5.4 -7 

Discussion The compliance of the components with drug-likeness (DL) 


In contrast to the conventional paradigm of "one drug, one 
target" in drug design, network pharmacology adopts a 
perspective of multi-targeted therapy to explore the intricate 
relationships between pharmaceuticals and diseases [17]. 
This unconventional approach leverages principles from 
systems biology, including network analysis, connectivity, 
and redundancy. Researchers have increasingly turned to 
network pharmacology studies [18-19] as a means to 
elucidate how medications interact with yet-to-be-identified 
signaling pathways. The Network Pharmacology framework 
[20] introduces novel insights into the systemic interplay 
between therapeutic targets and the entirety of a disease, 
rendering it a potent and promising tool for elucidating 
disease mechanisms at the systemic level and identifying 
potential bioactive ingredients. In the present study, a 
pioneering network was constructed to delve into the 
molecular processes associated with Nigella sativa. 


and oral bioavailability (OB) specifications was ascertained. 
Nigella sativa underwent further investigation, with 
emphasis on identifying bioactive targets related to 
colorectal cancer and their integration into a disease-centric 
network. The constructed network elucidated protein-protein 
interactions and various pathways, underscoring the 
potential of Nigella sativa bioactives in colorectal cancer 
inhibition. Key genes implicated in colorectal cancer, 
notably AKTI, IL6, ALB, HSP90AAI, ESRI, PTGS2, 
PPARG, MAPK3, GSK3B, and ERBB2, were identified 
through protein-protein interaction (PPI) analysis. 
Subsequent Gene Ontology (GO) and Kyoto Encyclopedia 
of Genes and Genomes (KEGG) analyses unveiled 
additional pathways, diseases, and disorders associated with 
these selected genes. Results from the GO enrichment study 
suggested a significant preventive effect of bioactives in 
colorectal cancer. Pathway analysis utilizing the KEGG 
database provided supportive evidence for the therapeutic 
potential of Nigella sativa in colorectal cancer treatment. IL- 


23 


Ansari Vikhar Danish Ahmad et al., JIPBS, Vol 11 (1): 10-24, 2024 


17 signaling, TNF signaling, Estrogen signaling, MAPK 
signaling, and Apoptosis emerged as significant nodes in the 
network. The findings indicated that Migella sativa 
influenced multiple signaling pathways, suggesting its 
applicability beyond prostate cancer, stomach cancer, and 
Th17 cell deficiency. Colorectal cancer, along with cancers 
of the colon, pancreas, endometrium, prostate, melanoma, 
bladder, lungs (small and large cell), liver, and stomach, 
demonstrated susceptibility to the anticancer effects of 
Nigella sativa. Validation of target accuracy was conducted 
through an additional docking port assessment. The 
evaluation of the components' affinity for their respective 
targets provided immediate insights into the structure- 
activity relationship. 


Conclusions 


This computational investigation utilizes network 
pharmacology and docking analysis to systematically 
explore the pharmacological mechanisms underlying the 
preventative and therapeutic effects of Nigella sativa in 
colorectal cancer. Emphasis is placed on elucidating the 
significant contributions of network pharmacology in 
unraveling the intricate mechanisms associated with Nigella 
sativa. The insights gained from this study may extend to 
benefit our understanding of analogous mechanisms in other 
related cancer types. As the current study has been 
conducted in line with the principle of network 
pharmacology but it is necessary to carry out validation 
study. 


Abbreviations: FAP, Familial adenomatous polyposis; 
HNPCC, Hereditary nonpolyposis colorectal cancer; MW, 
Molecular Weight; OB, Oral bioavailability; DL, Drug- 
likeness; BBB, Blood Brain Barrier; KEGG, Kyoto 
Encyclopedia of Genes and Genomes; GO, Gene Ontology; 
3D, Three Dimensional; PDB, Protein Data Bank; CAS, 
Chemical Abstracts Service; PPI, Protein Protein 
Interaction, DPED, | Duke's  Phytochemical and 
Ethnobotanical Databases; NS, Nigella sativa, CC, 
Colorectal Cancer; MD, Molecular Docking; SWP, Swiss 
Target Prediction. 
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